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Abstract. 

Measuring local velocities or entire flow rates in liquid metals or semiconductor 
melts is a notorious problem in many industrial applications, including metal casting 
and silicon crystal growth. We present a new variant of an old technique which relies 
on the continuous tracking of a flow-advected transient eddy current that is induced by 
a pulsed external magnetic field. This calibration-free method is validated by applying 
it to the velocity of a spinning disk made of aluminum. First tests at a rig with a flow 
of liquid GalnSn are also presented. 
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Compared to transparent flnids, liqnid metals and semicondnctor melts pose serions 
problems for the reliable determination of local velocities or integral flow characteristics 
[T]. Even the widely nsed nltrasonic Doppler velocimetry (UDV) faces difhcnlties when 
it comes to its application in very hot and/or chemically aggressive flnids, snch as liqnid 
steel or silicon. Fortnnately, the high electrical condnctivity that is responsible for the 
opaqneness of those flnids allows the ntilization of magnetic indnctive methods. 

These methods bear on applying magnetic helds to the flowing flnid and measnring 
appropriate featnres, e.g. amplitudes, phases, or forces, of the flow indnced magnetic 
helds. In the contactless indnctive how tomography (GIFT), entire three-dimensional 
how helds are reconstrncted from indnced held amplitndes that are measnred at many 
position aronnd the hnid, to which one or a few external magnetic helds are applied 
lai- A recently developed how rate sensor relies on the determination of magnetic 
phase shifts dne to the how mm- In the Lorentz force velocimetry (LFV) [5] one 
measnres the force acting on a permanent magnet close to the how, which resnlts as 
a direct conseqnence of Newton’s third law applied to the braking force acting by the 
magnet on the how. With this techniqne, it is now possible to measnre velocities of 
hnids with remarkably low condnctivities, snch as salt water [6]. 

A common drawback of all these methods is that the measnred signal is not only 
dependent on the songht how velocity, bnt also on the condnctivity of the hnid. Usnally, 
the signals are proportional to the so-called magnetic Reynolds nnmber Rm = fi^aVL, 
where po is the magnetic permeability constant, a the condnctivity of the liqnid, and 
V and L denote typical velocity and length scales of the relevant hnid volnme. In most 
cases, the signal depends also on geometric factors, so that the measnring system has 
to be calibrated anyway. Farther to this, the nse of permanent magnets (for LVF) or 
of magnetic yoke materials (for the phase-shift method) set limitations to the ambient 
temperatnre at the position of the respective sensors. 

The goal of the present paper is both to circnmvent the necessity to calibrate 
the measnrement system, and to mitigate the temperatnre limitation problem. The 
measnrement system to be presented belongs to a wider class of time-of hight methods 
that ntilize the existence of some traceable pattern in the hnid which is being advected 
by the how. By measnring the time of hight of the advected pattern between two 
positions along the how one can infer the how velocity. 

In one realization of this principle, the role of the pattern is played by some 
(nnspecihed) tnrbnlence elements moving close to the wall of the hnid. If the how 
direction is known, one can apply external magnetic helds at two positions and infer the 
how speed from correlations of the magnetic signals indnced by the tnrbnlence elements 
that are passing by [^. A similar correlation techniqne is presently nnder investigation 
in the time-of-hight LVF [8]. 

Another realization had already been developed in the 1960’s by Zheignr and 
Sermons j9j. In their ‘pnlse method’, the role of the traceable pattern is played by 
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an eddy current system that is imprinted into the conducting medium by switching on 
or off the current in one excitation coil, and by registering the flow advected pattern 
by another coil situated downstream the flow. The instant at which the current in this 
detection coil crosses zero indicates the passing-by of the center of the eddy current 
system. The flow velocity is then determined by dividing the distance between the two 
coils by the measured time interval. 

In |9], the method had been shown to allow to measure velocities of an aluminum bar 
from 20 m/s down to 1 m/s, but not signihcantly lower. One reason for this limitation 
was the exclusive focus on identifying the instant of vanishing current in the detection 
coil by means of an oscillograph. For small velocities (or, more exactly, for small Rm), 
when the decay of the eddy current is much faster then its advection, this instant smears 
out strongly so that its determination becomes less and less reliable. 

In the present paper, we will apply a modern data acquisition system and 
appropriate analysis methods in order to infer the velocity from the complete data sets 
of the induced current in three detection coils, rather than only one particular instant in 
one single coil. With view on measuring the entire transient signal after switching, we 
dub this method as ‘transient eddy current flow metering’ (TEC-FM). After explaining 
the measuring principle and technical realization, we will apply it to the determination 
of the velocity of a rotating aluminum disk, and to the flow of GalnSn in a circular tube. 
The paper closes with some outlook for future improvements and for possible industrial 
application. 

2. Measuring system and principle 

The heart of our measuring system consists of three small detection coils embedded 
into one larger excitation coil (see hgureH]). The excitation coil is installed parallel to 
the boundary of the fluid, so that the three detections coils are arranged in one row in 
direction of the flow. 

The block wiring diagram is shown in figure |2l The signal for the excitation coil is 
produced by an AGILENT 33200A waveform generator and then amplihed by a Rohrer 
PA2166 precision power ampliher. The induced signals measured by the three detection 
coils are amplihed by FEMTO DLPVA-IOO-B-D amplihers and then analyzed with a 
TEGTRONIX DP07104 oscillograph. 

In the following, we will focus on the case that the current in the excitation coil 
is on for t < 0, and is being switched off at t = 0 (note that the contrary case that 
the current is switched on at t = 0 is not completely equivalent: the subtle differences 
will be discussed in the outlook). With our excitation technique we achieve a smooth 
and non-oscillatory decay with a typical fall time of approximately 100 /rs. By this 
pulse, a ring-like eddy current system is induced in the nearby moving conductor. For 
a conducting half-space at rest, the evolution of this eddy current system can be given 
in a quasi-analytic form [10], which was also experimentally conhrmed in a liquid metal 
experiment m- Here, we illustrate the ring-like eddy current system for the case of a 
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Figure 1. TEC-FM sensor with cables (left), and detailed view on the coil system 
(right). The plugin connector (1) connects the power supply to the excitation coil (5) 
consisting of 300 windings of copper wire with 0.5 mm diameter on the base body (4). 
The voltages induced in the three detection coils (9) are distributed in the adapter (3) 
to three connectors (2). The detection coils consist of 35 windings with diameter 0.1 
mm that are wound on the winding former (8) that is fixed by screws (7) to the base 
body (4). The clamp (6) serves as a cable relief. 



Figure 2. Block diagram of the measurement system. 


solid electrical conductor moving to the right. Figure [3] shows the result of corresponding 
simulation with ANSYS multiphysics. The excitation coil is indicated below the body. 

This eddy current, in turn, produces a magnetic field b outside the fluid, whose 
63 ,-component (in flow direction) is zero exactly below the center of the (advected) eddy 
current ring. Close to this point, can be considered linear in x. If the detection 
coils are situated within this linear region, their signals can be exploited to identify the 
point of vanishing bx which represents the (moving) pole of the eddy current ring. The 
use of three coils enables to verify the linearity of the function bx{x). Yet, for the sake 
of simplicity, the following theoretical considerations will be restricted to the simpler 
two-coil system (hgurellj). 
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Figure 3. Simulated eddy current system induced in a right-moving solid body shortly 
after switching off the current in the emitter coil. The thick red arrows in the upper 
panel indicate the current lines, the thin gray lines are iso-contours of the current 
density. 



Figure 4. Schematic sketch of the evolution of the bx component of the magnetic field 
due to the eddy current. 


Shortly after switching off the current in the excitation coil, the component along 
the flow direction x can be parametrized in the following way (from here on, we replace 
bx by b): 

b{x,t) = m{t){x — Xp{t)) (1) 

with the decay function 


m{t) = toq exp {—t/ t) (2) 

where r is the typical decay time of the eddy current system in the liquid. Xp(t) 
represents the time-dependent pole position of the flow-advected eddy current system 
which can, in turn, be interpreted as the time integral over the velocity of the fluid: 

Xp{t) = [ Vx{t)dt . (3) 

Jo 

Hence, the velocity of the pole of the eddy-current related magnetic held can be 
determined as 

dxpit) 


Vp{t) := 


dt 


( 4 ) 


The pole position Xp at a given instant can be determined from the values b{xi,t) and 
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h{x 2 ,t) measured at the two detection coil positions Xi and X 2 , respectively: 


Xpit) = a:i - 


h{xifl){x2 - Xi) 
b{x 2 ,t) - b{xi,t) ■ 


(5) 


This system of equations would be sufficient if the magnetic helds were indeed measured 
by Hall or Fluxgate sensors. However, since we are using pick-up coils (also with view 
on later high-temperature applications), our measured signal is actually the voltage in 
the detection coil which is proportional to the time derivative of the magnetic helds 
rather than the magnetic held itself. 

The equation for this time-derivative of the magnetic held. 


b{x, t) = rh(t){x — Xp(t)) — m(t)xp(t) 


( 6 ) 


can be obtained from equation ([T]). The corresponding pole Xpp of this time-derivative 
is derived from setting the r.h.s. of equation ([6]) to zero, and utilizing the relations 
m{t)/m{t) = r, and Xp = Vxit)-. 


Xppit) = Xp{t) + TVxit) . 


(7) 


The velocity of this pole of the magnetic held derivative results than as 


Vpp := Xppit) = Vpit) + TVpit). ( 8 ) 

We see that the velocity of the pole of the derivative Vpp is identical to the velocity Vp 
of the pole of the held itself as long as the latter is time-independent. In our method, 
we will trace indeed the pole of the time-derivative, which is for every instant given by 


b{xxfl){x 2 - Xi) 

Xr,v[t) = Xl — - -^-, 

b{x 2 fl) - b{xifl) 

quite in analogy with equation ([5]). 


(9) 


3. Rotating aluminum disks 

In this section we will validate the TEC-FM method at a simple model for which the 
velocity is well known. For this purpose we chose a rotating aluminum disk of radius 
165 mm and thickness 10 mm, and apply our sensor at a radius of 145 mm (see hgure 
ED. 

The velocity at this radius can be pre-adjusted by the disk rotation rate. Figure [6] 
shows the time evolution of the signals at the three receiver coils, for 11 chosen velocities 
between 0 m/s and 5 m/s. 

From these data we estimate the position of the pole (hgure [7D for the hrst 2 ms 
after switching off the current. The black lines, drawn for comparison, would correspond 
to a perfect movement of the pole of the eddy current with the rotating disk. 

We see that, at least during the hrst 1 ms, the signal is very clean and allows to 
determine the speed of the pole. This is shown then in hgure [3 for which we have made 
an additional average over 10 pulsing events with an interval of 10 ms between them 
(so that we obtain one velocity value every 100 ms). Obviously, the velocity can be 
recovered with high accuracy. 
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Figure 5. Velocity measurement for a rotating aluminum disk (not to scale). 



Figure 6. Time evolution of the voltages (raw signals) measured in the three receiver 
coils (right) for 11 pre-adjusted disk velocities Vx (left). 


4. Flow of GalnSn in a pipe 

After having validated the method at a solid rotating disk, we apply it now to the case of 
a liqnid metal flow in a pipe. We consider a flow of the entectic alloy GalnSn (with the 
electrical condnctivity a = 3.3 x 10®) in a circnlar plastic pipe of inner radins 27.3 mm 
(for more details of the test rig, see [12]). Fignre [9]shows again the signal of the three 
coils for 15 different flow velocities (measnred by an independent flow meter). Note that 
the magnetic field decay is much faster than in the aluminum case with its much higher 
conductivity, since the decay time is roughly r = /iocrd^, with d a typical length scale 
of the order of the coil distance. This lower decay time would allow to provide velocity 
values approximately every 10 ms. 

Again, the time-dependence of the pole position Xpp (hgure 10) and the resulting 
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Figure 7. Time evolution of the reconstructed pole position Xpp of the eddy current 
for the 11 pre-adjusted disk velocities Vx- 



Figure 8. Inferred velocities Vpp of the aluminum disk for the 11 pre-adjusted velocities 
Vx- 


velocities VpP are shown (figure 11). Given the turbulent character of the pipe flow 
it is not surprising that the pole trajectory in figure 10 is much more noisy than the 
correponding pole trajectory for the aluminum disk. Consequently, we obtain also a 
larger deviation of the final velocity data (figure 11). Yet, the overall agreement with 
the pre-given velocity is quite reasonable, and might even be improved by a refined 
data analysis. 
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Figure 9. The same as Figure 6, but for the flow of GalnSn in a tube. 



Figure 10. The same as Figure 7, but for the flow of GalnSn in a tube. 


5. Outlook 

We have modified, and significantly enhanced, the pnlse method of Zheigur and Sermons 
[9] by nsing a compact design of three small detection coils embedded into a larger 
excitation coil, and by exploiting the complete signals measnred at these coils for the 
determination of the flow velocity. 

This ‘transient eddy cnrrent flow-metering’ (TEC-FM) has been validated at a 
rotating alnminnm disk. First experiments at a flow of GalnSn show the applicability 
of the method for liqnid metals. 

The method is calibration-free since it does not depend on the condnctivity of 
the moving body (at least, as long as the signals are still strong enongh). A fnrther 
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Figure 11. The same as Figure 8, but for the flow of GalnSn in a tube. 


advantage of the sensor system is the avoidance of any magnetic material. Indeed, only 
air coils are used which can easily be adapted even to very hot environments such as for 
slab casters or within the pullers of Czochralski crystal growth. 

A straightforward extension of the principle is to measure two dimensional flow 
structures close to the wall, just by adding two additional receiver coils orthogonal to 
the existing ones. 

As mentioned above, there are subtle differences between switching on and 
switching off the current. These have to do with the v x B induction term which adds 
a weak oo-shaped current system to the present ring-like one (as shown in hgure 3). A 
detailed study of these differences, and of their potential to infer also the conductivity 
of the fluid, is left for future work. 
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